Longitudinal and transverse spin asymmetries for inclusive jet production at 
mid-rapidity in polarized p + p collisions at y/s = 200 GeV 
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(The STAR Collaboration) 
We report STAR measurements of the longitudinal double-spin asymmetry All, the transverse 
single-spin asymmetry An, and the transverse double-spin asymmetries As and Att for inclusive jet 
production at mid-rapidity in polarized p + p collisions at a center-of-mass energy of y/s = 200 GeV. 
The data represent integrated luminosities of 7.6 pb -1 with longitudinal polarization and 1.8 pb _1 
with transverse polarization, with 50-55% beam polarization, and were recorded in 2005 and 2006. 
No evidence is found for the existence of statistically significant jet An, At,, or Att at mid- 
rapidity. Recent model calculations indicate the An results may provide new limits on the gluon 
Sivers distribution in the proton. The asymmetry All significantly improves the knowledge of gluon 
polarization in the nucleon. 
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I. INTRODUCTION 

Dccp-inelastic scattering (DIS) experiments with po- 
larized lepton beams and targets containing polarized 
nucleons have measured the inclusive spin structure func- 
tion gi(x,Q 2 ) of the nucleon over a wide range in 
Bjorken-x, 0.003 < x < 0.8, and Q 2 , 1 < Q 2 < 
100 GeV 2 /c 2 The DIS data, combined with data 

on the couplings in neutron and hyperon f3 decay, lead 
one to conclude that the quark contribution to the spin 
of a longitudinally polarized nucleon is only about 25%, 
well below naive expectations that the nucleon spin origi- 
nates mainly from the valence quarks. Perturbative QCD 
analyses [l8|]-[23j of the Q 2 dependence of gi(x,Q 2 ) gave 
the first insights into possible gluon polarization contri- 
butions, but the precision is thus far limited by the Q 2 
range that is accessible in the fixed-target experiments. 
Semi-inclusive DIS spin asymmetry measurements with 
identified pions and kaons have made it possible to de- 
lineate the quark and anti-quark spin contributions by 
flavor, and measurements with hadron pairs and open 
charm mesons have shown sensitivity to gluon polariza- 
tion d3-[2i|. 

Collisions of polarized proton beams at the Relativistic 
Heavy Ion Collider (RHIC) at Brookhavcn National Lab- 
oratory have made it possible to study proton spin struc- 
ture via hadroproduction of jets and other hard probes 
at ten-fold higher center-of-mass energies than previous 
DIS experiments. Of particular interest to the determina- 
tion of gluon polarization is the longitudinal double-spin 
asymmetry A LL , 



A 



LL 



(1) 



where a ++ and er H are the differential production cross 
sections when the beam protons have equal and opposite 
hclicities, respectively. STAR has published differential 
production cross section data for inclusive jet produc- 
tion at mid-rapidity with transverse momenta, px, in the 
range 5 < pr < 50 GeV/c [2!| that are well described by 
perturbative QCD calculations at next-to-leading order 
(NLO) [3(3. This supports the use of this framework in 
interpreting our measurements of All. 

Semi-inclusive deep-inelastic scattering experiments 
have also measured a broad ran ge of transverse spin 
asymmetries (see for example [3l| -|36j). including asym- 
metries sensitive to the Collins and Sivers effects. The 
Collins effect involves the convolution of quark transver- 
sity with the transverse-spin-dependent Collins frag- 
mentation function, which has been measured in e + e~ 
scattering [53, [Hj|. The Sivers effect ascribes a spin- 
dependent transverse momentum to the partons in a 
transversely polarized proton. Recently, global analy- 
ses have been performed to extract the quark transver- 
sity [1^] and parton Sivers [4(| distributions from the 
semi-inclusive DIS and e + e~ data. Measurements of 
the transverse double-spin asymmetry, Att, for inclu- 
sive jet production provide a complementary probe of 



quark transversity [4l|. It has also been proposed that 
the transverse single-spin asymmetry, An, for inclusive 
jet production may be sensitive to the Sivers effect [42j . 

In this article, we discuss the techniques that STAR 
uses to find and reconstruct inclusive jets in polarized 
p + p collisions, update our earlier analysis of All from 

2005 data [HI, and present precision data recorded in 

2006 on All. The All results significantly improve our 
knowledge of the gluon polarization in the proton. We 
also present the first results for the transverse single- 
spin asymmetry An and transverse double-spin asym- 
metries Ay, and Att (defined in Sect. IIV B| for the in- 
clusive production of mid-rapidity jets with transverse 
momenta up to 35 GcV/c in polarized proton collisions 
at y/s = 200 GeV. The An measurement may provide 
new limits on the gluon Sivers distribution in the proton 

m. 



II. EXPERIMENT AND DATA 
A. RHIC-STAR 

A schematic diagram of the STAR detector at RHIC 
is shown in Fig. [1] The detector studies collisions of 
independently polarized proton beam bunches, ranging 
from 50-112 bunches stored in each of two rings for a 
given fill. During the 2005 run, the proton bunches were 
loaded with alternating spin directions for the beam that 
circulated in the clockwise direction (blue) and with al- 
ternating spins for successive pairs of bunches for the 
beam that circulated in the counterclockwise direction 
(yellow). More complex 8-bunch polarization patterns 
were implemented in 2006 to further reduce possible sys- 
tematic errors associated with individual bunch patterns. 
Collisions at STAR are tagged by a coincidence of hits 
in the east and west Beam-Beam Counters (BBC) |44j . 
scintillation detectors that consist of 18 hexagonal tiles 
subtending the pseudorapidity interval, 3.4 < \rj\ < 5.0, 
on opposite sides of the interaction region. The perfor- 
mance of this detector as a luminosity monitor was cross- 
checked against a pair of hadronic Zero Degree Calorime- 
ters (ZDCs) located ±18 m from the detector center and 
at zero degrees relative to the beam axis. 

STAR employs several subsystems for particle track- 
ing and calorimetry (45l |. The central Time Projection 
Chamber (TPC) tracks charged particles with ^85% ef- 
ficiency over the full azimuth, < <fi < 2it, and pseu- 
dorapidity range, \rj\ <1.0, and falling to ^50% effi- 
ciency at \rj\ ~ 1.3. Surrounding the TPC is a Barrel 
Electromagnetic Calorimeter (BEMC) which consists of 
4800 Pb-scintillator towers covering the full azimuth over 
the range —0.98 < 77 < +0.98 in its final configuration. 
Each tower subtends an area of solid angle (A</> x A77) 
= (0.05 x 0.05). The positive r) side of STAR (west 
end) is covered by an additional 720 Pb-scintillator tow- 
ers comprising the Endcap Electromagnetic Calorimeter 
(EEMC), which extends calorimeter coverage over the 
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full azimuth for the pseudo-rapidity range 1.08 < 77 < 2.0. 
Both electromagnetic calorimeters are ~20 radiations 
lengths and ^1 strong interaction length deep. Fast sig- 
nals from the calorimeter towers are processed to classify 
triggers for events of interest. The reader is referred to 
Ref. for a comprehensive description of the STAR 
detector. 

The measurements presented here were taken over two 
different running periods during the years 2005 and 2006. 
In 2005, only the west half of the BEMC (0 < 77 < 0.98) 
was available. Between the data taking in 2005 and 2006, 
the BEMC commissioning was completed. This provided 
a more complete and robust picture of jets in our detector 
by doubling the acceptance and enabling measurements 
with jet cone radii, R = sj A<j) 2 + A77 2 = 0.7, which is 
larger than the value R = 0.4 that was used in our earlier 
analyses [H|4|. 



B. Triggers and Data Sets 

A minimum bias (MB) trigger was defined to be a co- 
incidence between any pair of BBC tiles from opposite 
sides of the collision region. This trigger has been shown 
to accept about 87% of the non-singly diffractive pp cross 
A redundant set of scalers recorded BBC 



44 



section 

tile and plane hits for each RHIC beam crossing, allow- 
ing the BBCs to be used as local luminosity monitors 
and polarimeters. The scaler system also recorded nu- 
merous combinations of hit conditions, including hits in 
the BBC on opposite sides of the interaction region in 
15 unequal intervals of the time difference between the 
hits. The intervals were chosen so as to give fine gran- 
ularity for beam-beam collisions that occurred near the 
center of the detector (z=0) and coarser granularity for 
events away from the center. We analyzed events from 



the intervals that correspond to a collision vertex selec- 
tion along the beam direction of approximately ±60 cm 
from the center of the detector. In this way, we matched 
the conditions for event selection with the conditions used 
in determining the relative luminosity for different spin 
combinations. The MB trigger was heavily prescaled to 
contribute only a few percent of the recorded data. 

Triggers for the selection of events with jets were con- 
structed by requiring substantial energy to be present in 
the BEMC. A High Tower (HT) trigger was defined by re- 
quiring a BBC coincidence plus at least one BEMC tower 
with a transverse energy greater than a given threshold. 
A Jet Patch (JP) trigger required a BBC coincidence, 
plus a transverse electromagnetic energy in a region of 
A77 x A<f> = 1.0 x 1.0 exceed a given threshold. 
The locations of the jet patches were fixed by hardware, 
with 12 such patches in the barrel calorimeter as shown 
in Fig. H 

In 2005, data were taken with a mixture of HT and JP 
triggers with different thresholds. The low (high) HT1(2) 
triggers required each accepted event to have at least one 
BEMC tower with transverse energy Et > 2.6(3.5) GeV. 
The low(high) JP1(2) trigger thresholds were set to Et > 
4.5(6.5) GcV. The HT1 and JP1 triggers were prescaled. 
There was considerable overlap among the triggers, with 
approximately half of the jets contained in the JP2 trig- 
ger sample. 

In 2006, the JP trigger was operated with a threshold 
of Et > 7.8 GeV early in the run, including the entire 
transverse polarization period. The threshold was then 
increased to 8.3 GeV for most of the longitudinal polar- 
ization period. Two additional triggers were also used. 
The first was a HT trigger that required Et > 5.4 GeV. 
The second was a refinement of the HT trigger (HTTP) 
that required a high tower to exceed a threshold of 3.8 
GeV, with an additional requirement of Et > 5.2 GeV 
in the 3x3 array of towers centered on the high tower. 
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The integrated luminosity for longitudinally polarized 
beam was 2.1 ptr 1 in 2005 and 5.5 ptr 1 in 2006. The in- 
tegrated luminosity for transversely polarized beam was 
1.8 pb" 1 during 2006. 



III. JETS 

Jets measure energy flow and are observable as com- 
posites of measured particle momenta. In the discus- 
sion below, we distinguish three categories: partem,, par- 
ticle, and detector jets. Our jet-finding and reconstruc- 
tion method uses the approach and algorithm adopted 
from Rcf. ;46j. Parton jets arc constructed from quarks 
and gluons either in theoretical calculations or in Monte 
Carlo simulated events prior to hadronization. Particle 
jets are constructed from the momenta of stable particles 
in Monte Carlo simulated events after hadronization. De- 
tector jets are constructed from real or simulated data. 
An important difference compared to many other exper- 
iments is that the charged hadron momenta are mea- 
sured in STAR with the TPC, rather than by a hadron 
calorimeter. 

We use comparisons of parton-to-particle and particle- 
to-detector jets to quantify the corrections needed to ac- 
count for contributions arising from spectator partons, 
effects of the underlying event in p + p collisions, unde- 
tected energy from, for example, neutrons, Kr, and neu- 
trinos, particle scattering out of the jet cone due to the 
hadronization process, bias and resolution effects intro- 
duced by our triggers and detectors, and uncertainties in 
the relative contributions of different partonic processes 
that result from uncertainties in the parton distribution 
functions. 



A. Jet-Finding and Reconstruction 

Detector jets are reconstructed in this analysis using 
a midpoint cone algorithm [4?} that combines charged 
tracks from the TPC and tower hits from the electromag- 
netic calorimeters. A primary vertex position is defined 
from the intersection of two or more charged tracks with 
the known transverse position of the beams. To be in- 
cluded in jet reconstruction, tracks are required to have 
a transverse momentum greater than 200 MeV/c, while 
tower hits must have a transverse energy Et exceeding 
200 MeV. Charged tracks are also required to contain 
>20 (>12) fit points in the TPC for 2005 (2006) and 
>51% of the maximum number of fit points allowed by 
the TPC geometry and active electronic channels. They 
are also required to have a distance-of-closest-approach 
to the primary vertex of dca < 3 cm. For the 2006 data, 
an additional px-dependent constraint was imposed on 
the transverse distance of the track from the beam line 
(dcaD): dcaD < 2 cm for p T < 0.5 GeV/c, dcaD < 
1 cm for pt > 1 GeV/c, and an interpolated cut in 
between. The tracking capabilities of the TPC gener- 



ally allow a determination of the dca value with a res- 
olution of 0.2-0.3 cm. These cuts are imposed to re- 
duce pileup and background tracks in the data. Par- 
ticles measured as TPC trajectories are assumed to be 
charged pions, whereas energy deposits in the electro- 
magnetic calorimeters are assumed to be photons. To re- 
duce double counting of energy, the average energy that 
a minimum-ionizing particle (MIP) would deposit is sub- 
tracted from the calorimeter tower energy if a TPC track 
points to the tower. If the calculated MIP energy depo- 
sition exceeds the energy observed in such a tower, the 
tower energy is discarded. 

The midpoint cone algorithm begins by collecting a 
list of tracks and tower hits with transverse momen- 
tum/energy greater than a set threshold value. These 
tracks and towers serve as seeds for the initial jet cones 
or proto-jets, defined as the collection of track and tower 
four-momenta inside of a cone of radius R, whose axis 
coincides with the i/^-weightcd centroid of the proto-jet 
four-momenta. Additional proto-jets arc formed from the 
midpoints between seeds and added to the list. At this 
point, a single track or tower may contribute to several 
proto-jets. 

Next, the algorithm decides whether to split or merge 
two proto-jets that have common four-momenta based on 
the fraction of energy shared by the two proto-jets. If the 
fraction is smaller than a specific value (0.5), the proto- 
jets are split into two jets and the shared four-momenta 
are assigned to the closest jets. If the fraction is greater 
than 0.5, the proto-jets are merged into a single jet. For 
2005 data, the jet cone radius was chosen to be R=0A. 
The cone radius was increased to R=0.7 for 2006 data 
to increase the efficiency and minimize the sensitivity to 
differences in the fragmentation of quark vs. gluon jets. 
Parameters used in the jet-finding algorithm are summa- 
rized in Table |IJ 



TABLE I: Midpoint Cone Algorithm Parameters 



Parameter 


2005 Jet-finding 


2006 Jet-finding 


Cone Radius (rad) 


0.4 


0.7 


Seed E T Threshold 


0.5 GeV 


0.5 GeV 


Assoc. Et Threshold 


0.1 GeV 


0.1 GeV 


Split/Merge Fraction 


0.5 


0.5 


Track p T Threshold 


0.2 GeV/c 


0.2 GeV/c 


Tower E T Threshold 


0.2 GeV 


0.2 GeV 


Jet rj 


0.2 < r\ < 0.8 


-0.7 < t] < 0.9 


Jet pt 


> 5.0 GeV/c 


> 5.0 GeV/c 



B. Event and Jet Cuts 

Events were removed from this analysis in the absence 
of a valid polarization measurement, relative luminosity 
value or BBC vertex information, or if the event failed the 
off-line verification of the trigger requirements. Events 
were also eliminated if the bunch identification tagged 
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them as originating from non-colliding, diagnostic (or 
"kicked") bunches. Kicked bunches are special bunches 
whose betatron orbits are deliberately amplified in order 
to give a large amplitude signal to a beam position mon- 
itor. Timing measurements of the kicked bunches give 
more precise measurements of the energy of the beam. 
After these requirements, longitudinal data samples of 
4.6M events were obtained in both 2005 and 2006. 

Many of the events contain two or more reconstructed 
jets. For the spin asymmetry measurements, only those 
triggered jets that contain a trigger tower or point to 
a triggered jet patch are included. Events that contain 
both a triggered jet and a non-triggered jet were included 
in the sample used to estimate the beam-gas background 
fraction (see next section). Further cuts, as described in 
the following sections, were imposed on the jets to elimi- 
nate backgrounds. After these requirements, the longitu- 
dinal data samples totaled 2.3(2. 1)M jets in 2005(2006). 
2% (4%) of the events contained two jets, both of which 
passed all cuts. 



C. Background Events 

1. Beam- Gas Events 

Energetic particles from beam-gas scattering and other 
non-collision background sources can pass through the 
beam-line shielding and then shower in the electromag- 
netic calorimeters. Most of these events lack signals in 
at least one of the BBCs, so they fail the trigger require- 
ment. Typically, they also fail our event cuts because 
they have no primary vertex within the active region of 
the TPC. However, such events can be mistaken for jets 
if they occur during beam crossings that also contain 
minimum-bias collisions. 

It is extremely rare to have two energetic background 
particles enter STAR concurrently, so non-collision back- 
ground events almost never contain two or more recon- 
structed jets. Also, background jets typically exhibit an 
abnormally large fraction of electromagnetic to total en- 
ergy, referred to as the neutral energy fraction (NEF), 
due to the lack of charged particle tracks that point back 
to an allowed vertex. We use these features to identify 
the non-collision background contribution to our jet sam- 
ple. 

Figure [3] shows the neutral energy fraction distribu- 
tions for our reconstructed jets when we divide them into 
two sub-samples. The di-jet events contain a second re- 
constructed jet with A<fi > ir/2. To enhance the di-jet 
statistics, the second jet is not required to satisfy the 
trigger independently. It is also allowed to fall within a 
larger pscudorapidity range than normal because it is not 
essential to reconstruct its energy precisely. The mono- 
jet events are the remainder. The measured neutral en- 
ergy distribution for di-jet events is well described by 
PYTHIA events processed through a GEANT model of 
the STAR detector (see Sect. MI D[) . The mono-jet events 
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NEF 

FIG. 3: (color online) Neutral energy fraction (NEF) distri- 
butions for mono- and di-jets for 2006 data (symbols) and 
Monte Carlo simulation (histograms). These events have un- 
corrected jet pt in the range 14.08 < pr < 17.31 GeV/c. 

are also well described except at large NEF, where a large 
enhancement due to non-collision background events is 
seen in the data. 




NEF 



FIG. 4: (color online) Fits of the mono/di-jet ratio vs. NEF 
for 2006 data. The total fit (red solid curve) includes a de- 
caying exponential for the background (blue dotted curve) 
and a first order polynomial for the signal (green dashed 
curve). These events have uncorrected jet pr in the range 
14.08 < p T < 17.31 GeV/c. 

The background component is determined from fits of 
the mono-jet to di-jet yield ratio as a function of NEF, 
as shown in Fig. 0] Signal events have a slowly vary- 
ing ratio, which is well described by a linear dependence 
in Monte Carlo simulations. The charged energy that is 
reconstructed in non-collision background events arises 
from tracks that accidentally point toward the calorime- 
ter energy deposition. These tracks arc typically associ- 
ated with the minimum-bias collision that was necessary 
to satisfy the trigger and produce a primary vertex. We 
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find that this component can be fit with an exponential 
function. 

The background mono-jets with NEF > 0.95 were used 
to set an upper limit on the asymmetry in the non- 
collision backgrounds. We then established the cut for 
signal events to minimize the quadrature sum of the sta- 
tistical and non-collision background uncertainties. For 
2006, we required NEF < 0.92. Less shielding was present 
in 2005, so we required NEF < 0.90 for all bins except the 
lowest px bin. It was difficult to isolate the background 
component in the fit for the lowest 2005 pt bin, leading 
to a very large uncertainty on any possible background 
asymmetry. Therefore, to be conservative we required 
NEF < 0.85 for this bin. 



2. Electron-like Events 

Initial studies of the NEF distributions for the HT and 
HTTP triggered events found an enhancement in the jet 
yield for NEF in the range 0.4 to 0.5 in both the data and 
Monte Carlo simulations. The enhancement, which was 
particularly prevalent at low jet pt, appeared primarily 
in the triggered mono-jet samples. The efficiency for ob- 
serving the second jet in a di-jct increased with jet pt, 
which indicated that the likely cause was events where 
the jet energy had been significantly overestimated. Fur- 
ther study determined that the enhancement arose from 
events where a conversion electron or positron that fired 
the HT or HTTP trigger also had a track reconstructed 
by the TPC. The jet finder will double-count the energy 
from such an electron or positron because it categorizes 
the TPC track and EMC tower as two high-energy par- 
ticles with similar 4- vector Et values. 

A set of cuts was implemented to minimize the recon- 
struction bias associated with the double counting. Jets 
in HT or HTTP triggered events were discarded if the 
highest pt track in the jet projected, within \Arj\ < 0.03 
and \A<j)\ < 0.027, to the location of the highest E T 
tower. The ratio of the tower energy to track momen- 
tum was required to be less than 1.2. Figure [3] shows 
that these cuts effectively removed these "electron-like" 
events. 



D. Simulations 

Monte Carlo events were generated using PYTHIA 
6.205 HI (2005 data) and PYTHIA 6.410 [3 [U 
(2006 data) with parameters adjusted to the CDF 'Tune 
A' settings [5l| and processed through the STAR detec- 
tor response package based on GEANT 3 [52| , including 
simulation of the trigger electronics. The PYTHIA pa- 
rameters corresponding to this tune are listed in Table [TT1 
In order to achieve a satisfactory simulation of the ob- 
served momentum balance for back-to-back jets (Fig. [5]), 
we incorporated an intrinsic parton transverse momen- 
tum of 1 GeV/c into the default PYTHIA model. Small 



discrepancies with the data for these model calculations 
may still be observed; however, for the purpose of esti- 
mating the systematic uncertainties, the shapes of the 
data distributions are reproduced sufficiently well by the 
simulations. We give further examples of the comparison 
of data and Monte Carlo in Figs. |6][8] for the jet pr spec- 
trum, track multiplicity and integrated transverse energy 
profile within the jets. These figures use 2006 data and 
Monte Carlo simulations. They require each jet to satisfy 
at least one of the HT, JP or HTTP trigger conditions, 
the same condition used in the asymmetry analysis. The 
effect of the triggers on the character of the jets is seen 
most readily in Fig. |H1 which plots the NEF distribu- 
tion of the jets for two different jet momentum ranges. 
At low momenta, the calorimeter-based triggers preferen- 
tially select jets with higher neutral energy fraction than 
jets at higher px- The systematic uncertainty associated 
with this trigger bias will be discussed below. 

TABLE II: CDF Tune A parameters. 



Parameter 


Value 


MSTP(51) 


7 


MSTP(81) 


1 


MSTP(82) 


4 


PARP(82) 


2.0 


PARP(83) 


0.5 


PARP(84) 


0.4 


PARP(85) 


0.9 


PARP(86) 


0.95 


PARP(89) 


1800 


PARP(90) 


0.25 


PARP(91) 


1.0 


PARP(67) 


4.0 



E. Jet Energy Scale 

Jet spin asymmetries are reported here as a function of 
the jet transverse momentum. However, a number of cor- 
rections must be made to the physical measurements to 
permit comparisons to the parton-level cross sections and 
jet pt definitions used in theoretical calculations. Some 
corrections are best subsumed into the systematic uncer- 
tainties on the asymmetries themselves, while others are 
more naturally applied as shifts to the jet momenta. 

1. TPC and Calorimeter Calibration 

The TPC calibration proceeds in several steps j53j . 
First, the drift velocity, which is monitored during run- 
time via laser ionization of the gas, is determined to 
an accuracy of approximately 0.03%. Then, distortion 
of the tracks due to misalignments of the readout sec- 
tors, drift distortions in the magnetic field and construc- 
tion imperfections are removed by examining extremely 
rigid tracks in different volumes of the TPC. Finally, 
the effects of space charge due to positive ion build- 
up at high event rates is monitored and corrected by 
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FIG. 5: (color online) Relative difference in pr for back- 
to-back di-jets for 2006 data. Both jets satisfy all cuts, 
including the trigger-matching requirement. These events 
have uncorrected jet pr in the range 14.08 < (pt,i + 
Pt.2)I1 < 17.31 GeV/c. For this case, we find ~ 5% dis- 
crepancy between data and Monte Carlo simulations. Essen- 
tially all the bins match simulations to within 10%. The px 
resolution is estimated to be ~ 23% from this graph. 




30 35 40 45 50 
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FIG. 6: (color online) Raw jet yield versus uncorrected jet 
transverse momentum in 2006 data (points) compared with 
Monte Carlo simulations (histogram). 



examining the distance of closest approach to the pre- 
sumed common vertex of an ensemble of tracks. This 
latter quantity is monitored function of instanta- 
neous luminosity throughout the runs, and the correc- 
tions are updated periodically. The resulting hit er- 
rors on the tracks are in the range 300-550 fim. The 
overall momentum resolution of the TPC tracks is ap- 
proximately A Pt /pt ~ 0.01 + 0.005p T /(GcV/c) for 
Pt < 10 GeV/c [S3|- The TPC tracking efficiency is 
~ 85 ± 5%. 

The electromagnetic calorimeter was calibrated using 
a sample of identified electrons in the TPC data with 
momenta between 1.5 and 15 GeV/c that satisfied strict 
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FIG. 7: (color online) Top Panel: Average track multiplicity 
in the reconstructed jet versus jet pr for 2006 data (solid 
symbols) and Monte Carlo simulations (open symbols). Ver- 
tical bars represent the r.m.s. width of the multiplicity dis- 
tributions rather than the uncertainty on the mean. Bottom 
panel: Ratio of data to Monte Carlo simulations. Vertical 
bars show the statistical uncertainties. 
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FIG. 8: (color online) Fraction of the total jet transverse 
energy found within a cone of radius AR centered on the 
reconstructed thrust axis, illustrating the jet profile. Data 
from 2006 (symbols) and Monte Carlo simulations (curves) 
are shown for three different corrected jet pr bins. 



geometrical and isolation conditions. The extrapolated 
trajectory was required to remain completely within a 
tower and the maximum energy of the nearest neighbors 
was limited. The tower calibration was then determined 
from the distribution of the ratio of energy observed in 
the towers divided by the momentum of the track. The 
variation of this ratio as a function of distance from the 
center of the tower was studied in data and Monte Carlo, 
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FIG. 9: (color online) Neutral energy fraction (NEF) for 2006 
data (symbols) and Monte Carlo simulations (histograms) for 
uncorrected jet pt in the ranges 7.6 < pr < 9.3 GeV/c (upper 
panel) and 39.6 < pr < 48.7 GeV/c (lower panel). These 
plots demonstrate the substantial bias introduced at low pr 
by triggering on only electromagnetic energy. 



tector jets, bin-by-bin and for each trigger type, that are 
then combined to produce the pt shift for the data. The 
main systematic uncertainties on this px shift are due to 
assumptions about the proportions of different partonic 
subprocesses contributing to the jet spectrum. The un- 
certainties were estimated by recalculating the px shifts 
for different subprocesses in PYTHIA, taking the max- 
imum deviation for any of them, then adding the sta- 
tistical uncertainty in quadrature with this value. The 
smaller cone size used in 2005, compared to 2006, results 
in a larger asymmetry in the associated systematic error. 
Table EH (2005) and Table [TV] (2006) give these correc- 
tions for both years. For each transverse momentum bin 
(first column) we give the mean detector jet pr (second 
column) and the corrected mean particle jet pr (third 
column). The fourth column in each table gives the es- 
timated uncertainties on the pt shift. The fifth column 
arises from uncertainties in jet energy scale due to pos- 
sible inaccuracies in the calibration and performance of 
the TPC and EMC detectors. The comparison of particle 
and detector jets was repeated for a variety of calibration 
ranges, tracking inefficiencies, and detector states in or- 
der to estimate conservatively the range of possible effects 
on the pt shift due to detector performance. 



and a correction was applied to compensate for this vari- 
ation. The correction amounts to a maximum of 8% at 
the edge of the fiducial cut. 

The total uncertainty in the calorimeter response to 
jets is estimated to be 4.8% for the present analysis. 
This includes the uncertainty from the electron calibra- 
tion plus contributions from the uncertainty in the dif- 
ferential response of the calorimeter to hadronic vs. elec- 
tromagnetic energy and the ability of the Monte Carlo 
simulation to describe the precise light output from the 
scintillators when energy is deposited very close to the 
edge between two towers. The nominal full scale gain for 
the individual calorimeter towers was changed from 28 
GeV in 2005 to 60 GcV in 2006. Therefore, the calorime- 
ter gain uncertainties are independent between the two 
years, even though the fractional uncertainties are equal. 



2. Jet pt Scale Corrections 

Our largest pt correction accounts for the difference 
in particle and detector jet pt scales. The combination 
of a steeply-falling pt dependence in the jet yield and 
a jet transverse momentum resolution of ~ 23% (seen in 
Fig. [5J causes substantial bin migration. Thus, on av- 
erage, lower momentum jets are reconstructed as higher 
momentum jets. The most straight-forward method of 
dealing with this effect is to apply a pt shift to correct 
the average value of the detector jet pt within a bin. 

This correction is calculated by comparing the pt cen- 
troid for Monte Carlo simulations of particle versus de- 



3. Pile-up Corrections 

Pile-up refers to the rate-dependent correction for 
charged tracks and calorimeter hits that were acciden- 
tally added to a jet during reconstruction. The largest 
pile-up contribution came from out-of-time tracks that 
were nonetheless reconstructed within the 40 /is TPC 
read-out period. Additional sources included multiple 
events within the same bunch crossing or beam halo back- 
ground that was coincident with a hard collision. 

To estimate the size of this correction, during nor- 
mal data-taking a small fraction of events were taken 
with a random trigger, that is, data taken during nom- 
inal beam crossings, but with no detector requirements. 
These events are expected to contain the effects of the 
pile-up energy alone, including the correct averaging over 
the instantaneous luminosity during the data-taking. 

Jets were reconstructed in a sample of normal events. 
The tracks and calorimeter hits from random events in 
the same run were then added to these same normal 
events and a new set of jets were reconstructed. Separate 
average shifts were calculated for tracks and calorimeter 
hits and these average shifts were applied to the final 
jet pt spectra. For 2005(2006) this amounted to a shift 
of 0.008(0.050) GeV/c per jet. The larger correction for 
2006 reflects the larger cone size used in the jet recon- 
struction and the higher instantaneous luminosity that 
was available. 
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TABLE III: Jet transverse momentum bins and corrections for Y2005 data. For each pr bin, the average detector jet pr and 
the corrected particle jet pr values are listed. Also shown are the evaluated uncertainties on these pr values arising from the 
PT-shift procedure, from the residual uncertainties in our detector simulations, and from the uncertainties on our hadronization 
and underlying event estimation and the unknown QCD scale. The sum in quadrature of these uncertainties is tabulated in 
the final column. All uncertainties are in units of GeV/c. 



Measured pr Range 


< Pt > 


Corrected < pr > 


p T Shift 


Detector Sim. 


Hadr/UE/QCD 


Total 


(GeV/c) 


(GeV/c) 


(GeV/c) 


Uncertainty 


Uncertainty 


Scale Uncertainty 


Uncertainty 


5.00 - 6.15 


5.58 


5.32 


+0.18/-0.24 


±0.23 


±0.27/-0.15 


±0.40/-0.36 


6.15 - 7.56 


6.86 


6.30 


+0.11/-0.09 


±0.25 


±0.37/-0.28 


±0.46/-0.39 


7.56 - 9.30 


8.43 


7.06 


+0.36/-0.04 


±0.27 


±0.43/-0.34 


±0.62/-0.44 


9.30 - 11.44 


10.37 


8.67 


+0.20/-0.16 


±0.35 


±0.52/-0.39 


±0.66/-0.55 


11.44 - 14.08 


12.76 


10.73 


+0.09/-0.07 


±0.41 


±0.68/-0.54 


±0.80/-0.68 


14.08 - 17.31 


15.70 


13.08 


+0.08/-0.07 


±0.52 


±0.75/-0.54 


±0.92/-0.75 


17.31 - 21.30 


19.31 


16.00 


+0.19/-0.22 


±0.63 


±0.80/-0.50 


+1.04/-0.83 


21.30 - 26.19 


23.75 


19.39 


+0.30/-0.33 


±0.77 


±0.94/-0.53 


+1.25/-0.99 


26.19 - 32.22 


29.21 


23.57 


+0.38/-0.29 


±0.94 


+1.12/-0.61 


+1.51/-1.16 


32.22 - 39.63 


35.92 


28.07 


+0.58/-0.40 


±1.12 


+1.29/-0.65 


±1.80/-1.36 



TABLE IV: Jet transverse momentum bins and corrections for Y2006 data. Details are the same as given in the caption for 
Table [m] 



Measured pr Range 


< PT > 


Corrected < pr > 


pr Shift 


Detector Sim. 


Hadr/UE/QCD 


Total 


(GeV/c) 


(GeV/c) 


(GeV/c) 


Uncertainty 


Uncertainty 


Scale Uncertainty 


Uncertainty 


7.56 - 9.30 


8.43 


8.51 


±0.52/-0.37 


±0.35 


±0.51 


±0.81/-0.72 


9.30 - 11.44 


10.37 


10.32 


±0.30/-0.35 


±0.40 


±0.57 


±0.76/-0.78 


11.44 - 14.08 


12.76 


12.17 


±0.25/-0.23 


±0.46 


±0.76 


±0.92/-0.92 


14.08 - 17.31 


15.70 


14.41 


±0.08/-0.08 


±0.55 


±0.81 


±0.98/-0.98 


17.31 - 21.30 


19.31 


17.15 


+0.21/-0.17 


±0.66 


±0.86 


+1.10/-1.10 


21.30 - 26.19 


23.75 


20.45 


±0.13/-0.16 


±0.80 


±0.98 


+1.27/-1.28 


26.19 - 32.22 


29.21 


24.42 


±0.10/-0.12 


±0.97 


±1.17 


+1.52/-1.52 


32.22 - 39.63 


35.92 


29.41 


+0.22/-0.27 


±1.17 


±1.37 


+1.82/-1.82 


39.63 - 48.74 


44.19 


34.72 


±0.90/-1.22 


±1.38 


±2.03 


+2.61/-2.74 



4- Hadronization and Underlying Event Corrections 

The pt shift applied above is still not complete. There 
is a further difference in scale between parton jet mo- 
mentum and particle jet momentum that may be divided 
into two partially compensating effects: underlying event 
(UE) and out-of-cone (OOC)[55( fragmentation. 

The OOC effect causes a reduction of the measured 
Pt in particle jets due to fragmentation and hadroniza- 
tion of the parent parton outside of the jet cone. This 
correction may be sub-process dependent as quarks are 
expected to have a harder fragmentation spectrum than 
gluons. The underlying event causes an increase in the 
measured px in particle jets due to the inclusion of parti- 
cles arising from interactions between spectator partons 
in the proton remnants and from additional hard par- 
tonic scatterings in the event. The UE is expected to be 
isotropic in r\ x <f> space and has been found to be largely 
independent of jet pr- 

The combined effect from the OOC and UE on the jet 
Pt scale was estimated by comparing reconstructed jets 
at the fragmented parton (FP) and particle (PART) level 
in events generated by the PYTHIA 6.3 Monte Carlo 
package O [5(| with parameters set to the CDF Tunc 
A values. Jets at the FP stage contain only the frag- 



mented partons resulting from the scattered partons and 
the initial and final radiation (MSTP(61), MSTP(71)). 
At the FP stage the underlying event and hadroniza- 
tion (MSTP(81), MSTP(lll)) are turned off. Note that 
MSTP(81) only controls the multiple parton interaction 
component of the UE and does not include effects from 
remnant interactions. Jets at the PART level contain the 
stable, hadronized, final-state particles resulting from the 
interaction in addition to any initial and final-state radi- 
ation. The reconstructed jet pt scale at the PART level 
in simulations is comparable to the experimentally mea- 
sured jet scale after corrections for detector resolution 
and trigger bias are included. 

The total change in jet pt scale, Apt — Pt P —Pt ART : 
depends on the radius of the jet cone. Generally the shift 
is smaller for jets reconstructed with the larger cone ra- 
dius in 2006. The reduced shifts at larger pt for R=0.7 
indicates that OOC effects become less important as the 
cone radius increases. The lower pt behavior is domi- 
nated by UE effects. As expected, the Apt was found 
to be sub-process dependent and larger for gluon jets. 
Therefore these effects have been included as a system- 
atic uncertainty, instead of a correction, on the measured 
detector±trigger corrected jet pt- These uncertainties 
are given m the sixth columns of Table Hill and Table HVl 
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IV. SPIN ASYMMETRY ANALYSIS 

A. The Spin Asymmetry All 

Experimentally, the double longitudinal spin asymme- 
try defined in Equation 1 was evaluated according to: 

E(PiP2)(N++-rN+-) 

All = o , 

J2(P 1 P 2 ) 2 (N++ + rN+-)' 

where Pi 2 denote the measured beam polarizations, 
N ++ and denote the inclusive jet yields for equal 

and opposite proton beam hclicity configurations, respec- 
tively, and r is the ratio of measured luminosities for the 
two helicity configurations. Each sum is performed over 
runs that last from 10 to 30 minutes so that the mea- 
surements are sampled on time-scales faster than typical 
variations in the beam polarizations and relative lumi- 
nosities. 



1. Beam Polarization 

The beams are injected and circulated as bunches in 
the RHIC rings with their spins oriented in the verti- 
cal direction. Their polarizations are measured using 
Coulomb-nuclear interference polarimeters [57l . |58| that 
are calibrated against a polarized gas jet target [59( lo- 
cated at other interaction regions around the RHIC ring. 
The magnitude of the polarization is measured and mon- 
itored throughout the beam stores from these locations 
and is generally in the range from 50-55%, with a statis- 
tical uncertainty of ~ 1%. 

2. Relative Luminosity 

The relative luminosity for each polarization combina- 
tion in a run was calculated from the sum of BBC coin- 
cidences over a run, after sorting bunches for each spin 
combination. Since these rates enter directly into the 
expression for the asymmetries, care was taken to ensure 
these data were consistent and systematically understood 
to a level commensurate with the size of the asymmetry 
being measured. 

The BBC and ZDC analog pulses are discriminated 
and the coincidence signals are used as input signals to a 
time-to-amplitude converter, whose output is converted 
to a 4 bit time difference signal. For each beam cross- 
ing (every 106.5 ns), 4 bits for the BBC coincidence sig- 
nal, 4 bits for the ZDC coincidence, and 7 bits for the 
beam crossing number are distributed to a set of 4 re- 
dundant scaler boards. These scaler data are then ex- 
amined for statistical consistency; in general we find ex- 
cellent agreement among all methods of luminosity mea- 
surement. However, a small fraction of runs (< 1%) were 
rejected due to inconsistencies among the BBC measure- 
ments from different electronics channels. 



As a further safeguard against detector failure or sub- 
tle physical effects, the BBC relative luminosity mea- 
surements were cross-checked against the ZDC measure- 
ments. The BBCs are sensitive to the total non-singly 
diffractive pp cross section by intercepting single charged 
particles over a broad rapidity range at moderate trans- 
verse momenta [44| . This hypothesis has been supported 
by PYTHIA calculations used in the design of these 
detectors, and by direct measurements (Vernier scans), 
which demonstrate that a very large fraction of the non- 
singly diffractive cross section is indeed measured [60l ]. 
The ZDCs detect neutral particles like neutrons and 7r°s 
close to the rapidity of the beam, and are thus sensitive 
to types of collisions which arc very different from those 
sampled by the BBCs. The small acceptance of the ZDCs 
limits the statistical precision of this comparison (for the 
present data); however, this has the additional advantage 
of allowing an examination of rate-dependence of the lu- 
minosity measurements as well. The result of this com- 
parison was consistent between the two years and gives 
a conservative systematic uncertainty on the relative lu- 
minosity of slightly less than 10~ 3 . 

Corrections to the luminosity are expected due to ac- 
cidental coincidences and under-counting of multiple in- 
teractions in a single beam crossing, as explained in 
Ref. [6l|. These effects have been examined for the 
relative luminosities encountered in the 2005/2006 runs 
and found to be negligible compared to the systematic 
uncertainty assigned to the relative luminosity by com- 
paring BBC and ZDC measurements. 

B. Transverse Spin Asymmetries 

In the 2006 run, STAR recorded 1.8 pb" 1 of jet 
data from transversely polarized proton-proton collisions. 
These data have been analyzed in order to measure the 
transverse spin asymmetries A^, A-£, and Att- 

In a coordinate system where the positive z axis and 
pseudorapidity are defined by the momentum direction 
of beam 1, the polarization directions of the beams are 
along the y axis, and the azimuthal angle <j) is defined 
relative to the x axis, we can write the jet production 
cross section for the two transversely polarized protons 
as in Ref. (62|: 

= 1 + PiP 2 • A^(rj,p T ) 
+ cos(0) ■ [Pi • A N (n,p T ) - P 2 • A N (-r),p T )) (3) 
+ PiP 2 -cos(20) • A tt (v,Pt) 

Additional transverse spin asymmetries can be defined 
related to particle correlations within a jet [U [63| that 
are beyond the scope of this paper. 

The extraction of the asymmetries proceeded as fol- 
lows. The single-spin asymmetry, An, was determined 
by combining the spin directions for one beam to approx- 
imate an unpolarized 'target'. The single-spin asymme- 
try for each beam was determined separately, using the 
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FIG. 10: (color online) Ajv as a function of the corrected mean 
Pt for 2006 transverse data. The panels present An for four 
different 77 bins. An is the left-right single-spin asymmetry for 
a transversely polarized beam. The errors shown combine the 
statistical uncertainties, which dominate, with all systematic 
uncertainties except trigger and reconstruction bias. See Sect. 
IV Al for a discussion of the latter. 



cross-ratio technique [64[, and the results combined. The 
statistical precision for this measurement was sufficient 
to allow measurements as a function of jet transverse 
momentum in 4 bins of pseudorapidity relative to the 
polarized beam, as shown in Fig. 1101 

The double-spin asymmetry Ay, was determined by av- 
eraging over the entire range of pseudorapidity and az- 
imuth, using analysis procedures identical to those for 
All. The results for Ay, are given in Fig. [TT] We are not 
aware of any theoretical predictions for Ay ■ However, it 
plays an important role in the estimation of the system- 
atic uncertainty in All due to residual transverse spin 
components in the beam (see Sect. IIV C 21) . A statis- 
tically significant measurement of Att could be made 
only by averaging over the full data set and extract- 
ing the coefficient of the cos(20) dependence on the az- 
imuthal angle from the fit shown in Fig. [TJ] We find 
A T t = -0.0049 ± 0.0046. This precision is not yet 




FIG. 11: (color online) At, versus corrected mean px for 2006 
transverse data, averaged over the range \r/\ < 0.8. As, the 
transverse double-spin asymmetry for transversely polarized 
beams, is defined in more detail in the text. 
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FIG. 12: (color online) Measured transverse double-spin 
asymmetry versus azimuthal angle cj> for 2006 transverse data. 
The points are fit with the function As +Att cos(2<^>). These 
data are averaged over 7.5 < pr <40 GeV/c and \q\ < 0.8. 



sufficient to confront predictions of Att due to quark 
transversity |4lj . 



C. Systematic Uncertainties 

1. Trigger and Reconstruction Bias 

Jet pt resolution effects shown in Fig. [5] cause aver- 
aged shifts in the jet pt scale as discussed in the previ- 
ous section. Additionally, jet events arc selected based 
on neutral energy triggers, that preferentially select jets 
with characteristics that differ from those of the unbiased 
jet distribution. For example, for the same jet momen- 
tum, the HT trigger will preferentially fire on jets with 
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FIG. 13: (color online) Trigger and reconstruction bias esti- 
mates for the longitudinal double-spin asymmetry All versus 
jet pr for a representative set of polarized parton distribution 
models [I},[6g]-[62||. 

a high-energy leading particle while the JP trigger will 
fire on jets with larger radii. The relative proportions of 
quark-quark, quark-gluon and gluon-gluon interactions 
in an unbiased sample of events are fixed predominantly 
by the (well-determined) unpolarizcd structure functions. 
However, effects of the trigger and jet reconstruction can 
change the relative proportions in our measured sample, 
and this can bias our measurements of A LL and A N . 

Calculation of the biases introduced by our trigger and 
jet reconstruction further depends on assumptions of the 
polarized parton structure functions. Our calculations 
must account for the uncertainty in polarized gluon con- 
tributions and, to a lesser degree, the associated uncer- 
tainties in polarized quark and sea contributions. 

Parameterizations of the polarized parton distribution 
functions are combined with PYTHIA parton kinematic 
variables to generate predictions of All vs. px specific 
to a particular model at both the particle and detector 
levels. A broad range of polarized parton distribution 
functions were adopted for these calculations. Eventu- 
ally, only those that predict distributions for All vs. pt 
consistent with the general trends of our measured results 
were included in the bias estimate. 

The method of calculating the trigger and reconstruc- 
tion bias is then as follows. Relative fractions of jets 
resulting from HT, JP, and HTTP triggered events have 
been measured in data and found to agree with simula- 
tions. The logical OR of the three trigger types is there- 
fore used to measure the trigger bias. For each px bin 
and polarized parton distribution model, All is calcu- 
lated both at the particle and detector levels. Detector 
All points are shifted as discussed above to correct to 
the particle jet pr scale. The residual difference between 
the particle jet All and the shifted detector jet All 
represents the bias for that model and pr bin. To be 
conservative, the systematic uncertainty for each bin is 



assigned to be the largest positive and negative differ- 
ence of all the allowed models. In Fig. [13] we show the 
result of this calculation for five representative models. 
The GRSV +0.3 and GRSV -0.3 models (HIH use thc 
GRSV functional form for the gluon polarization, with 
the integral fixed at the two respective values. DSSV 
[13] is a recent fit that includes input data from RHIC, 
in addition to DIS and semi-inclusive DIS. GS-C J20j is 
an early model that has a large gluon polarization at low 
x. These five models span a range in gluon polarization 
that is wider than permitted by our results. The sys- 
tematic uncertainties as a function of pt are listed in the 
second column of Tables [V] and IVII 



2. False Asymmetries from Residual Transverse Spin 
Effects 

False asymmetries that mimic our All signal can arise 
from a combination of physical and experimental sources. 
To obtain longitudinal collisions at STAR, the trans- 
versely polarized beams are rotated to the longitudinal 
direction, then back again to transverse on either side of 
the interaction region by a pair of helical dipoles known 
as spin rotators. Inaccuracies in the adjustment of the 
spin rotator currents leave small transverse components 
for both beams in the collision region. Our transverse 
spin asymmetry measurements allow us to put stringent 
limits on the associated false asymmetries. 

The transverse asymmetries in the central rapidity re- 
gion are expected (and measured) to be small. How- 
ever, in the presence of non-longitudinal polarization 
components, the asymmetry A^ can contribute directly 
to the observed All signal. Local measurements of the 
transverse polarization components of both beams during 
longitudinal running were made by examining thc single- 
spin asymmetries observed in BBC tile hits. The trans- 
verse single-spin asymmetry (An) has been reported pre- 
viously for this detector [3 [HI • As given in these ref- 
erences, it is in the range of ~ 6-7 xl0~ 3 and can be 
calibrated to a high accuracy during transverse running. 
Because the BBC is a highly segmented detector, com- 
binations of up/down and left/right scatterings can be 
used to measure the transverse polarization components 
for both beams. These measurements were made contin- 
uously through the data-taking. The residual transverse 
components for both beams during the nominally longi- 
tudinal run were weighted by integrated luminosity for 
different periods of adjustment for the spin rotators. De- 
noting the angle of the polarization with the longitudinal 
axis as 9, values of tan((9) between 0.02-0.18 were mea- 
sured, with an average magnitude equal to ~ 0.1 for both 
beams. 

Because measurements of As were consistent with 
zero, we do not make a correction to All for this con- 
tribution, but instead assign a systematic uncertainty 
on our All measurements. We combined the measure- 
ments of the transverse polarization components with the 
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TABLE V: pr-dependent systematic uncertainties for Y2005 data. The trigger and jet reconstruction bias, non-longitudinal 
beam polarization, and beam-gas background systematic uncertainties on the measured 2005 All are given. 



Pt 


Trigger bias and Jet Recon. 


Non-longitudinal Pol. 


Beam Gas Background 


(GeV/c) 


(X 10" 3 ) 


(x 10" 4 ) 


(x KT 4 ) 


5.32 


-1.67/+1.67 


±4.97 


±13.50 


6.30 


-1.47/+1.29 


±2.49 


±8.07 


7.06 


-1.44/+1.10 


±1.54 


±7.97 


8.67 


-1.80/+2.66 


±1.20 


±8.40 


10.73 


-1.24/+2.32 


±1.16 


±6.82 


13.08 


-1.34/+2.50 


±1.31 


±5.34 


16.00 


-1.88/+2.86 


±1.74 


±4.55 


19.39 


-2.74/+2.74 


±2.69 


±3.51 


23.57 


-3.91/+3.91 


±4.75 


±0.50 


28.07 


-3.44/+5.60 


±9.56 


±0.00 



TABLE VI: pr-dependent systematic uncertainties for Y2006 data. Details are the same as given in the caption for Table FVl 



Pt 


Trigger bias and Jet Recon. 


Non-longitudinal Pol. 


Beam Gas Background 


(GeV/c) 


(x 10" 3 ) 


(x 10" 4 ) 


(x 10" 4 ) 


8.51 


-2.00/±3.37 


±0.90 


±15.17 


10.32 


-1.07/±1.97 


±0.71 


±7.66 


12.17 


-1.26/±1.99 


±0.68 


±4.92 


14.41 


-0.58/+1.11 


±0.77 


±3.43 


17.15 


-0.43/±0.70 


±1.02 


±3.57 


20.45 


-0.72/+1.52 


±1.58 


±4.48 


24.42 


-1.03/±3.92 


±2.79 


±7.52 


29.41 


-1.57/+5.46 


±5.61 


±7.92 


34.72 


-2.88/±6.93 


±12.85 


±5.01 



4- Polarization and Relative Luminosity 

Systematic errors arising from beam polarization and 
relative luminosity measurements are treated separately 
due to their correlated effects on the data. An error in the 
relative luminosity measurement would result in a shift 
of the All data points by an additive constant, while an 
error in the measurement of the polarization magnitude 
would scale the magnitude of the All data. Therefore 
we quote these quantities separately. 

The systematic uncertainty on the relative luminos- 
ity was determined by comparing the BBC and ZDC 
measurements, which were found to be consistent at the 
< 1(T 3 level. The BBC-ZDC difference was used to esti- 
mate the possible size of the uncertainty on All due to 
the errors on the relative luminosity as 8 All ~9x 10 -4 
for both 2005 and 2006 data. 

The fractional systematic uncertainty for the quan- 
tity P1P2 quoted by the RHIC CNI polarimeter group 
(common to all RHIC experiments) is 9.4%(8.3%) for 
2005(2006) [6^. Polarization measurements from differ- 
ent years have contributions which may be identified as 
cither uncorrclatcd or correlated. The total error was 
conservatively estimated by assuming the latter portion 
to be 100% correlated from year to year. The correlated 
error in the normalization of beam polarizations comes 
mostly from an unpolarized molecular hydrogen back- 



uncertainties on the measurements of Ay, in each mo- 
mentum bin to give a conservative (maximal) estimate 
of the systematic uncertainty: 

SA L l = |tan(0i) tan(0 2 ) cos(0i - ^2) x A s \ 

~ \6i6 2 5Av\. (4) 

where ^1(^2) and <f>i{4>2) are the polar and azimuthal an- 
gles of the polarization vectors for beams 1 (2), respec- 
tively. In keeping with the spirit of estimating this un- 
certainty conservatively, the value of cos{4>\ — ^2) was 
set equal to unity. The uncertainty on All due to non- 
longitudinal components of the beams as a function of 
Pt is listed in the third column of Tables [V] and IVII 



3. Beam Gas Background 

The systematic uncertainty on the residual beam back- 
ground was conservatively estimated to be the larger of 
either the measured effect or the statistical uncertainty of 
the measured effect on the asymmetry. These numerical 
values range from ~ 15 x 10~ 4 in the smallest pt bins to 
less than half this value at higher px- This uncertainty as 
a function of px is listed in the fourth column of Tables IVl 
and ED 
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ground in the gas jet polarization measurement. These 
errors represent an overall scale uncertainty on our mea- 
surements, and are common with the polarization uncer- 
tainties of concurrent measurements performed by the 
PHENIX experiment. They are therefore quoted sepa- 
rately to facilitate comparison of different data sets and 
to identify correlated errors, where possible. 



D. False Asymmetries 

All measurements were examined as a function of time 
to ensure the absence of non-statistical variations. In 
addition, the data for different beam and bunch com- 
binations were combined to form parity-violating single 
and double spin asymmetries. These are expected to 
be highly suppressed and provide internally consistent 
cross-checks on the validity of the measurements. Dou- 
ble spin asymmetries were formed from the "like-sign" 
and "unlike-sign" bunch combinations, and single spin 
longitudinal asymmetries were formed for each beam di- 
rection. No false asymmetry was found to be significantly 
different from zero. 



V. RESULTS AND DISCUSSION 



A. An Results 



independent evaluations of the asymmetries and system- 
atic uncertainties. While the list of corrections from both 
years are the same, individual items differ in magnitude 
and range from year to year. Furthermore, the different 
jet definitions and pseudorapidity intervals for the two 
years lead to different expectations for All from model 
calculations. We therefore do not combine the results 
from different years, but present them separately. Our 
results from year 2005 data are given in Fig. 03] and Ta- 
ble |Vnl The corresponding results for year 2006 data are 
given m Fig.rjS]and Table EEED Note the different scales 
on the vertical axes for the two figures. 
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In a recent model calculation, D'Alesio et al. [42j con- 
clude that An for inclusive jets at mid-rapidity arises 
solely from the gluon Sivers effect. Within their model, 
the leading contribution to the systematic uncertainty 
on our measured inclusive jet An arises from the fact 
that our triggers have different efficiencies for detecting 
jets from cither quark-quark, quark-gluon, or gluon-gluon 
scattering. Although the latter two processes dominate 
at low jet pt, the detector efficiency is largest for quark- 
quark scattering. This can lead to a measured An that 
is smaller in magnitude than the true value. We have 
used our Monte Carlo simulations to estimate the size of 
this effect. We find that correcting for this bias would 
increase the magnitude of An by up to 40% for \ow-pr 
jets, dropping to 25% at 15 GeV/c and 15% at 30 GcV/c. 

D'Alesio et al. find that the current upper limit on the 
gluon Sivers distribution would lead to \An\ of 4-5% at 
Pt = 8 GeV/c, dropping to -2.5% at 15 GeV/c Q. 
The results in Fig. [TU] indicate that our measured An 
is substantially smaller than these upper limits. Thus, 
they may provide new constraints on the magnitude of 
the gluon Sivers distribution in the proton. 

B. All Results 

The different detector geometries, triggers, jet defini- 
tions, and measured pseudorapidity ranges of the two 
different data taking years reported here have demanded 



FIG. 14: (color online) All for inclusive jet production versus 
corrected jet px for 2005 data. For 2005 data we used a cone 
radius equal to 0.4 and a pseudorapidity range for the jet 
thrust axis of 0.2 < r\ < 0.8. The error bars are statistical. 
The gray boxes show the systematic uncertainties. 



TABLE VII: The final measured All and pr values from the 
2005 data sample. Data cover the range 0.2 < 77 < 0.8 with 
a jet cone radius of 7?=0.4. Statistical and systematic uncer- 
tainties are listed for All (note the systematic uncertainty is 
asymmetric as described in the text). The p-r-dependent cor- 
rections listed in Table IVl were combined in quadrature with 
the systematic uncertainty of 9 x 10~ 4 in relative luminosity 
to give the total. The pr values shown are the results after 
applying all pr corrections discussed in the text. 



Pt 


A LL 


stat. err. 


sys. err. 


(GeV/c) 


(x 10" 3 ) 


(x 10" 3 ) 


(x 10" 3 ) 


5.3 +0.4/-0.4 


5.3 


±5.9 


+2.4/-2.4 


6.3 +0.5/-0.4 


-2.7 


±5.4 


+1.8/-1.9 


7.1 +0.6/-0.4 


2.4 


±5.7 


+1.7/-1.9 


8.7 +0.7/-0.6 


14.3 


±6.7 


+2.Q/-2.2 


10.7 +0.8/-0.7 


-6.7 


±8.7 


+2.6/-1.7 


13.1 +0.9/-0.8 


2.6 


±12.7 


+2.7/-1.7 


16.0 +1.0/-0.8 


-14.6 


±20.3 


+3.0/-2.2 


19.4 +1.3/-1.0 


-52.2 


±35.0 


+2.9/-2.9 


23.6 +1.5/-1.2 


56.9 


±67.1 


+4.1/-4.1 


28.1 +1.8/-1.4 


146 


±138 


+5.7/-3.7 
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FIG. 15: (color online) All for inclusive jet production versus 
corrected jet pr for 2006 data. For 2006 data we used a 
jet cone radius equal to 0.7 and a pseudorapidity range of 
— 0.7 < < 0.9 to reflect the increased detector coverage. The 
error bars are statistical. The gray boxes show the systematic 
uncertainties. 



TABLE VIII: The final measured All and pr values from the 
2006 data sample. Data cover the range —0.7 < r\ < 0.9 with 
a jet cone radius of R—0.7. Details are the same as given in 
the caption for Table IVTTI 



PT 


All 


stat. err. 


sys. err. 


(GeV/c) 


(x 10" 3 ) 


(x 10" 3 ) 


(x 10" 3 ) 


8.5 +0.8/-0.7 


2.7 


±5.3 


±3.8/-2.6 


10.3 ±0.8 


3.3 


±4.3 


+2.3/-1.6 


12.2 ±0.9 


9.9 


±4.1 


+2.3/-1.6 


14.4 ±1.0 


1.2 


±4.5 


+1.5/-1.2 


17.2 ±1.1 


9.2 


±5.8 


+1.2/-1.1 


20.5 ±1.3 


25.7 


±8.8 


+1.8/-1.3 


24.4 ±1.5 


25.6 


±15.4 


+4.1/-1.6 


29.4 ±1.8 


22.0 


±31.0 


+5.6/-2.1 


34.7 ±2.6/-2.7 


12.0 


±70.6 


+7.1/-3.3 



C. Comparison to Theory 

The theoretical curves shown in Figs. [14] and [15] arc de- 
rived from NLO calculations of spin asymmetries based 
on the code of Jager et al. [30|. This code provides both 
the polarized and unpolarized proton-proton cross sec- 
tions for an input cone of radius R centered at rapidity 
y and averaged over azimuth to 0(3) in a s . Expressions 
for all 2 — > 2 and 2 — > 3 processes were derived analyti- 
cally in a small cone approximation, with the subsequent 
integrals evaluated numerically using a Monte Carlo ap- 
proximation. The results were compared to more com- 
plete calculations [70[ also using Monte Carlo evaluation 
of the integral, but without the small cone assumption 

The code allows as input the colliding energy of the 
protons, the jet cone radius, the jet rapidity and the 
jet transverse momentum intervals for integration of the 
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FIG. 16: (color online) Scale dependence of the relative asym- 
metry (defined in the text) as a function of jet transverse mo- 
mentum. The asymmetry is calculated for the GRSV Stan- 
dard set of parameters. The black (red) lines are the ratios 
of the asymmetries calculated for half (twice) the nominal 
renormalization and initial state scales of (j,r — — pr- 
The dotted curves are for a jet cone radius of 0.4 and solid 
and dashed curves are for a cone radius of 0.7. 



cross section. The code requires assumptions for two 
scale inputs: the initial-state factorization scale (/xj) 
and the renormalization scale (fin). We take the value 
A'/ = M-R = Pt of the jets. The polarized and unpolar- 
ized inclusive jet cross sections are calculated by separate 
programs. 

As originally configured, the programs require polar- 
ized and unpolarized parton distribution functions that 
are sampled by the Monte Carlo portions of the pro- 
gram over a wide range of x and Q 2 . These are tab- 
ulated at fixed values beforehand and interpolated to 
the required precision. The original configuration also 
allowed for a selection among several sets of unpolar- 
ized (CTEQ5, CTEQ6M, CTEQ6M.1) and polarized 
(GRSV2000 STND, MAX, MIN, ZERO) parton distri- 
bution functions. All calculations here use the pdf set 
CTEQ6M [7ll and a cone radius of 0.4 or 0.7, as noted. 

The possible size of higher order corrections to the 
cross sections is conventionally estimated by varying the 
factorization and renormalization scales by a factor of 
2 about the nominal scale. In Fig. [TH1 we plot calcu- 
lations that show the scale-dependence of the relative 
asymmetry, A l l/ A LL , as a function of jet transverse 
momentum for both jet cone radii. In this expression, 
A° LL is the asymmetry calculated with the nominal scales 
Hi. r = Pt, while All is calculated for scales of 2px and 
Pt/2. The parton distribution functions used in this cal- 
culation were the CTEQ6 set and GRSV2000 STND. In 
general, the larger the cone radius, the less sensitive is 
the calculation of the spin asymmetry (not cross section) 
with respect to scale and higher order contributions. The 
choice of cone radius equal to 0.4 for 2005 was made in 
consideration of the size of the (instrumented) part of the 
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detector, in order that acceptance uncertainties would 
not dominate the systematic uncertainty on the jet en- 
ergy. With increased EMC coverage in year 2006, our 
sensitivity to scale variations is lessened. 

The impact of these data on previous constraints of 
the integral of the polarized gluon distribution function, 
AG, from deep-inelastic lepton scattering data is evident 
in Figs. ITU and ITS! The measured All values are seen 
to lie in the region below the previous best-fit DIS curve, 
GRSV standard [3(1 [sS HH ■ The remaining curves asso- 
ciated with GRSV are the polarized parton distributions 
re-fit to constrain AG to a series of values spanning the 
full range —g(x) < Ag(x,Qo) < g(x), that is the gluon 
spins may be fully polarized in either direction, or in- 
terpolated to intermediate values using a common func- 
tional form. 

Data from a single experiment cover a limited kine- 
matic range in x and Q 2 , making the measurement of 
the total integral AG(Q 2 ) at a specific Q 2 impossible. A 
rigorous extraction of AG requires the incorporation of 
these inclusive jet asymmetries, along with other RHIC, 
DIS and SIDIS polarized scattering data, into a global 
analysis. For example, the AAC analysis [72| demon- 
strated that while the inclusion of PHENIX pion lon- 
gitudinal double-spin asymmetries [73| available at that 
time had only a small influence on the optimum fit, the 
uncertainties on the gluon polarized parton distribution 
function were significantly reduced over the fit obtained 
using DIS data alone. 

This type of analysis, more recently performed by 
de Florian et al. (DSSV) [H uses NLO pQCD fits 
to the world data set (including the STAR 2005 [H and 
a preliminary version of the 2006 jet asymmetries pre- 
sented here), constrained by a functional form that de- 
fines Ag(x, Q 2 ) in the unmeasured regions of x space, to 
extract the spin-dependent parton densities. The DSSV 
global analysis is based on Mellin moments, which allows 
a certain amount of x integrated data to be included in 
the fits. This is an especially important development 
for RHIC data, where the statistical precision within our 
kinematical constraints thus far only allows the examina- 
tion of AG over a limited range of 0.02 < x < 0.2. The 
DSSV best fit finds the gluon polarization to be much 
smaller than that in GRSV standard throughout the x 



region which is currently constrained by data. Further- 
more, in the x region sampled by RHIC data the DSSV 
X 2 + 2% upper limit on Ag(x,Q 2 ) at Q 2 = 10 GcV 2 
is roughly half the GRSV standard value (see Fig. 2 in 



VI. SUMMARY 

In summary, we have reported an analysis of spin de- 
pendencies in the inclusive production of mid-rapidity 
jets with transverse momenta up to 35 GcV/c in polarized 
p+p collisions at y 7 ! = 200 GcV from data recorded in 
2005 and 2006. No evidence is found for the existence of 
statistically significant transverse asymmetries Ajq, A-%, 
and Att- The An result may provide new limits on the 
gluon Sivers distribution in the proton. The longitudi- 
nal double-spin asymmetry All has been compared with 
NLO perturbativc QCD evaluations based on selected po- 
larized parton distribution functions to demonstrate its 
sensitivity to the value of the gluon helicity distribution 
inside the proton. 
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